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Abstract 

The effective field theory approach is applied to the three-nucleon process of S = 
1/2 neutron-deuteron scattering in the S-wave, including the effective range parameters 
summed at all orders. This is achieved through a modification of the power counting, 
following a recent suggestion in the literature. It is shown that with this procedure, 
while the convergence of the loop integrals improves, one cannot meaningfully include 
a three-body effective term to describe low-energy data affected by the presence of the 
triton bound state. 



I. Introduction 



Recently, three-body interactions in nucleon systems have attracted considerable attention, 
both within more classical approaches based on potential models [|T[ and in the context of 
effective field theories |, 1 (EFT's) §, |, |, |7|, |. 

The three-nucleon system is a natural test-ground for the understanding of nuclear forces 
that has been reached in the two-nucleon system. A convenient example is given by neutron- 
deuteron scattering H || ||, since in this case Coulomb interactions are negligible. Neutron- 
deuteron scattering involves two S- wave channels, corresponding to total spin 5 = 3/2 and 
5 = 1/2. In the 5 = 3/2 channel all spins are aligned and the two-nucleon interactions are 
only in the 3 Si partial wave. The two-body interaction is attractive but the Pauli exclusion 
principle forbids the three nucleons to be at the same point in space. Therefore this channel 
is insensitive to short-distance physics and disallows a three-body bound state. Thus one is 
able to obtain precise (~ 4%) predictions in a straightforward way [§, ||, ||. In the 5 = 1/2 
channel two-nucleon interaction can take place either in the 3 Si or in the 1 So partial waves. 
This leads to an attractive interaction which implies a three-body bound state, the triton. 
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The S = 1/2 channel also shows a strong sensitivity to short-distance physics as the Pauli 
principle does not apply. 

In this note EFT is applied to the latter channel, but including the effective range parame- 
ters for the 3 Si and 1 Sq NN scattering partial waves in the calculation, following a suggestion 



proposed in |H|. In that paper the authors argued that taking both the scattering length 



a and the effective range r in the NN scattering effective range expansion of order 
and summing range corrections to all orders, improves convergence and may solve problems 
encountered in the pionful theory [D] with KSW power counting [[|. The authors then suc- 
cessfully applied this modified power counting to a few two-body processes. It is therefore 
important to test this new scheme on a three-body system. It will be shown that the new 
power counting improves the convergence of the loop integrals necessary for the nd scattering 
calculation, so that the cutoff introduced in [|, § is no longer needed for the system of inte- 
gral equations without a three-body interaction term. However, even this improved ultraviolet 
behavior of the kernels in the integral equations is not enough, since adding a three-body in- 
teraction introduces new divergences which make the equations insensitive to the three-body 
term in the Lagrangian, as it will be shown below. 

This independence of the amplitudes with respect to an additional parameter prevents 
us from applying the fitting procedure of ref. || and therefore from predicting the energy 
dependence of the nd phase shifts. 



II. Formalism 



It is convenient to use a Lagrangian JT^, |9| expressed in terms of two auxiliary fields d % and 
t A with the quantum numbers of the deuteron and of a dibaryon field in the x So channel of 
NN scattering respectively: 
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where N = ( p ) is the nucleon doublet of two-component spinors, the subscripts d and t denote 
the 3 Si and 1 Sq channel of NN scattering, and M is the average nucleon mass. P l d and P A 



are the projectors onto the 3 Si and 1 S channels, respectively 



(2) 



with cr (r) the Pauli matrices acting in spin (isospin) space. The coefficients yt, A\ ^ and 
A^ encode all short distance physics - like pion and u exchanges, quarks and gluons, and 
resonances like the A. Here it is necessary to split both A's into leading (A^ -1 )) and subleading 
pieces (A (0) ). 

Since the theory is nonrelativistic, all particles propagate forward in time, the nucleon 
tadpoles vanish, and the propagator for the nucleon fields is 



iS(p) 



p -p 2 /2M + ie 



(3) 
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The deuteron and dibaryon propagators are more complicated because of the coupling to two- 
nucleon states. For instance, the bare deuteron propagator is simply a constant, — i/A^ 1 , but 
the full propagator gets dressed by nucleon loops to all orders as illustrated in Fig. [I]. 



Figure 1: The deuteron propagator at LO from the Lagrangian ([J). The thick solid line 
denotes the bare propagator , the double line its dressed counterpart. 

The nucleon-loop integral has a linear ultraviolet (UV) divergence which can be absorbed in 
yJ/A^ 1 , and a finite piece determined by the unitarity cut. Summing the resulting geometric 
series leads to the deuteron (and analogously the dibaryon) propagators in terms of physical 
quantities at LO: 

.l~i fiv 

iA?(p) = i^A d (p) = 



MVd i d -^/ p --M Po -ie 
4vri 5 AB 
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7d = \JMB = 45.7066 MeV and the deuteron binding energy is B = 2.225 MeV. The typical 
momentum j t = 1/df of the virtual bound state is extracted from the scattering length in this 
channel, a t = —23.714 fm. 

III. S = 1/2 nd Scattering 

The Lagrangian (|I|) can now be used to describe the nd-scattering in the S = 1/2 channel. The 
treatment involves coupled channel equations both in the 3 Si and 1 So partial waves. The spin 
zero dibaryon field t also contributes in intermediate states of nd amplitudes. It is possible to 
obtain a system of two coupled integral equations (previously derived using a different method 
flT3"| ) for the d + N — > d + N amplitude it^(k,p, e) and for the d + N ^ t + N amplitude 
if t A (k,p,e), pictorially represented in Fig. 0. 

A momentum cutoff A has to be introduced in the integral equations. In the limit A — > oo 
these equations do not have a unique solution because the phase of the asymptotic solution 



is undetermined ||14j| . For a finite A this phase is fixed and the solution is unique. However, 
the equations with a cutoff display a strong cutoff dependence that does not appear in any 
order in perturbation theory. The amplitude td(p,k = const.) shows a strongly oscillating 
behavior ||. This cutoff dependence is not created by divergent Feynman diagrams. It is 
a nonperturbative effect and appears although all individual diagrams are superficially UV 
finite. For an extended discussion on this effect, see 0. 

The solution is to add one-parameter three-body force counterterm H(A)/A 2 that runs 
with the cutoff A M. This counterterm represents a three-body force which is obtained by 
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Figure 2: The coupled set of Faddeev equations which need to be solved for t& t &t LO in the 
doublet channel. The double dashed line denotes the dibaryon Geld t\ A . 
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in the Lagrangian (|T|). This yields the equations: 
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71 



dq q 2 t d (k, q) 



K(p, q) + 
K(p, q) + 



2H(A) 

2H(A) 
3A 2 



o 

A 



¥ - ME - is - ld 



7T Vt 



dq q 2 t t (k,q) 



2f - ME - is - lt 



(6) 



t t {k,p) 



l2y d y t M 



K(p,k) 



x 



x 



K(p, q) + 
K(p, q) 4 



2H{A)~ 

2H(A) 
3A 2 



2H{A) 
3A 2 



3 yt 
nyd 



7l 



dq q 2 t t (k,q) 



*£-ME-ie- lt 



dq q 2 t d (k,q) 



*3--ME-ie- ld 



(7) 



where k (p) denote the incoming (outgoing) momenta in the center-of-mass frame, ME 
3k 2 /4 — 7J is the total energy, the kernel K(p, q) is given by 



1 ( q 2 +pq — p 2 — ME 
K [p, q) = - — In 



2pq \q 2 — pq — p 2 — ME J ' 
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and H(A) is defined as follows: 



H(A) 



sin[s ln(A/ A*) - arctg(l/ s )} 
sin[s ln(A/A*) + arctg(l/s )] 



(9) 



In eq. (f|) sq rs 1.0064 and A* is a dimensionful parameter that determines the asymptotic 
phase of the off-shell amplitude |4j], and is fitted to reproduce the experimental value for the 



S = 1/2 nd scattering length, 4 1/2) = (0.65 ± 0.04) fm pL yielding A, = 0.9 fm 



IV. Effective Range Expansion 

The cutoff dependence is expected to be eliminated from the equations if the effective range 
parameters p d and rot are included in the effective range expansion and are treated like the 
scattering length a, i.e. taken of order Q~ l , and summed to all orders. Ref. |]10) gave reasons 



to conclude that this power counting should improve the overall convergence. In the 3 Si 
channel, the expansion around the deuteron pole gives 

fc cot (5) = - ld + ^p d (k 2 + 1 2 ) + ... (10) 

for the nucleon-nucleon phase shifts. In the singlet channel 1 So, no real bound state exists, so 
the condition to impose is that the effective range expansion 

fccot(£) = -- + K 0t k 2 + ... (11) 

is satisfied. Resumming the effective ranges pd and rot in the full deuteron and dibaryon 
propagators generates 
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The numerical values p d = 1.765 fm and r ot = 2.73 fm [T(| have been usedQ. 

It is still possible to introduce in the computation a three-body term, which now is inde- 
pendent from the cutoff A and is in fact a constant contact interaction. The resulting integral 
equations are completely analogous to (H), (0), but now with the propagators ( ]T2|) . The inte- 
gration is carried on the full real positive axis without any cutoff. The integral equations are 
then solved numerically using the techniques outlined in ||. 

The solutions for the amplitudes coming from the equations do not show any visible depen- 
dence on the three-body interaction term H . It is possible to extract the energy dependence 
for A; cot (5), where k is the incoming momentum of the particles in the center of mass frame, 
from the behavior of the on-shell amplitude t d (k, k). Results are plotted on Fig. |^. Compare 
with the results given in ref. || (Fig. p. 

^or an early computation including effective ranges and comparisons with potential models, see p7| . 
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Figure 3: Energy dependence for 5 = 1/2 
nd scattering with H = or with H ^ 
and A — > oo when the effective ranges are 
included in the calculation. 
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Figure 4: Cutoff-insensitive energy de- 
pendence for 5 = 1/2 nd scattering for 
A* = 0.9 fm -1 , as given in ref. ||, with- 
out resumming the effective ranges. Data 
are from the phase-shift analysis of van 
Oers and Seagrave [18|] (crosses) and a 
measurement by Dilg et al. |15[] (dia- 
mond) . 



The explanation for this outcome stems from the fact that, while the convergence of the 
equations has improved because of the inclusion of the effective ranges (which multiply a factor 
oc q 2 in the denominator of the propagators), this is not yet sufficient to make the equations 
completely convergent. After subtracting from the amplitudes td,t(k,p) the parts satisfying 
the (convergent) equations without the three-body term H, the remaining parts t' dt (k,p) are 
described by equations containing linear divergences (plus convergent terms). Regularizing 
these divergences by introducing a cutoff A drives t' dt (k,p) to approximate quantities oc 1/A, 
i.e. to negligibly small numbers if the cutoff is set sufficiently high. td,t{k,p) are then effectively 
cutoff- independent, but independent from H as well. 



V. Conclusions 

The study of three-nucleon systems using EFT methods in the 5 = 1/2 channel is more com- 
plicated than for the 5 = 3/2 channel. While for the latter nd scattering accurate predictions 
are obtained ||, the 5=1/2 channel displays a strong cutoff dependence even though all 
individual diagrams are UV finite. This dependence can be eliminated at LO only for the 
equations without three-body interactions, if the effective range parameters, obtained from 
the effective range expansion in AW scattering, are taken into account in the deuteron and 
dibaryon propagators. This is achieved modifying the usual power-counting scheme, assuming 
pd and rot of order Q~ x like the scattering length a. 

Unfortunately this improvement is not sufficient to eliminate all the infinities originated in 
the loop diagrams from the integral equations. The remaining divergent terms have the effect 
of driving the three-body-dependent parts of the amplitudes to negligible values, while the 
amplitude parts dependent only on the finite equations (those without the three-body term) 
prevail. 

This makes eqs. (0), (H), after the inclusion of the effective range parameters, insensitive 
to the three-body term necessary to describe the influence of the triton bound state and to 



reproduce the experimental data for nd scattering in the S-wave doublet channel. 

Alternatively, one can argue that the three-body force term arise not at LO but at NLO. 
Yet in this case it must generate a contribution to the amplitude much bigger than the ~ 30 % 
variation one expects going from LO to NLO. 

Therefore this test of the new power-counting procedure does not yield the successful 
results achieved in the two-body problem case. Yet it is possible that a better treatment of 
this problem implies the explicit inclusion of the triton propagator in the theory, analogously 
to what has been done so far for the deuteron and the dibaryon, and according to the spirit 



ofref. 10 
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